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Abstract 

We describe the nature of charge transport at non-zero temperatures (T) 
above the two-dimensional (cQ superfluid-insulator quantum critical point. We 
argue that the transport is characterized by inelastic collisions among ther- 
mally excited carriers at a rate of order ksT/h. This implies that the trans- 
port at frequencies u <C ksT/fi is in the hydrodynamic, collision-dominated 
(or 'incoherent') regime, while lo S> ksT/h is the collisionless (or 'phase- 
coherent') regime. The conductivity is argued to be e 2 /h times a non-trivial 
universal scaling function of hoo/ksT, and not independent of Tiw/A^T, as has 
been previously claimed, or implicitly assumed. The experimentally measured 
d.c. conductivity is the hydrodynamic hu/ksT — > limit of this function, 
and is a universal number times e 2 /h, even though the transport is incoher- 
ent. Previous work determined the conductivity by incorrectly assuming it 
was also equal to the collisionless ftw/fceT — ► oo limit of the scaling function, 
which actually describes phase-coherent transport with a conductivity given 
by a different universal number times e 2 /h. We provide the first computa- 
tion of the universal d.c. conductivity in a disorder-free boson model, along 
with explicit crossover functions, using a quantum Boltzmann equation and 
an expansion in e = 3 — d. The case of spin transport near quantum critical 
points in antiferromagnets is also discussed. Similar ideas should apply to 
the transitions in quantum Hall systems and to metal- insulator transitions. 
We suggest experimental tests of our picture and speculate on a new route to 
self-duality at two-dimensional quantum critical points. 
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I. INTRODUCTION 



The charge transport properties of systems near a zero temperature quantum phase tran- 
sition have been the subject of a large number of experimental studies in the last few years. 
Systems in two spatial dimensions have been of special interest, and important examples 
are the superfluid-insulator transition in disordered thin films [0-0 and Josephson junction 
arrays ||, the transitions from quantized Hall plateaus a quantum critical point in 



the doped cuprate compounds 0, and a recent quantum transition in Si MOSFETs JTO 
In three dimensions, the metal-insulator transition in doped semiconductors has seen many 
years of study (JTI] , and evidence for scaling collapse of data near the quantum critical point 



has finally begun to emerge ||12|| . We shall present the results of this paper using the lan- 
guage of the superfluid-insulator transition in thin films in zero external magnetic field, but 
we believe our ideas are much more general and also have significant implications for the 
other systems noted above. We shall also discuss the extension of our results to the case of 



spin transport near quantum critical points in two-dimensional antiferromagnets [13| . 

A particular focus of the experiments on the superfluid-insulator transition has been the 
theoretical prediction, made in the seminal work of Fisher et al. |14|] and Cha et al. ]T5[ , that 
the conductivity at the critical point in two dimensions is non-zero and equals the quantum 
unit of conductance {e 2 /h) times a universal number. However there remains an impor- 
tant dichotomy between experimental and theoretical studies of the superfluid-insulator 
transition in disordered thin films, which is crucial to our discussion. All experiments are 
performed at a low, but non-zero, temperature T and have measured (d.c.) conductivities 
at frequencies, u, which easily satisfy hu -C k B T. In contrast, all of the theoretical analyti- 



cal |T4H20|| and exact diagonalization |2T[ work has computed the universal conductivity at 
T = 0, which is equivalent to the regime hu ^> k B T . Numerical Monte Carlo work has also 
measured the universal conductivity [ p!5| , p2} -p5| , but it involves an analytical continuation 
from imaginary Matsubara frequencies i27rpksT ' jh (p > 1, integer) to real frequencies: such 
a procedure is insensitive to any structure in the conductivity at hu <C k B T, and so these 
measurements are effectively also in the regime hu 3> k B T. This discrepancy, between the- 



ory and experiment, in the orders of limit of u — > and T — > was noted by Cha et al. [T3 



and Wallin et al. ||25|| , but they asserted that the conductivity approaches a universal value, 
at low temperatures and frequencies, which is independent of the value of the ratio hu/k B T. 

In this paper, we shall argue that the critical conductivity is in fact not independent of 
hu/k B T, and is instead a highly non-trivial, but universal, function of hu/k B T; this shall be 
explicitly demonstrated in a computation of the crossover function in a model of disorder- 
free bosons. We shall show that the physics of the hu 3> k B T and hu <C k B T regimes are 
quite distinct, and shall provide, for the first time, a model computation of the universal 
properties in the limit hu <C k B T relevant to d.c. transport measurements. We hope that 
our results will stimulate experimental measurements of the a.c. conductivity at frequencies 
which can explore the crossover at hu ~ k B T '. 

Our arguments rely on the physical picture of the T > dynamics at quantum-critical 
points developed in the context of quantum antiferromagnets p6|JT3|j28|j29|] . It was argued 
in these works that the order parameter dynamics is relaxational, with a relaxation rate 
~ k B T/h. Using the orientation of the order parameter to define a phase, we then obtain 
a phase relaxation rate l/r^ ~ k B T/h, and this rate is analogous to the phase relaxation 
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rate of disordered metals (the latter rate is much smaller than k B T '/% in the disordered 
metal [30], but becomes of order k B T/h at the metal insulator transition-see Section |VD|). 



At time scales of order r v , the relaxational dynamics cannot be described by effective classical 
models of the types discussed in Ref |31|, as a necessary condition for classical behavior is 



that the relaxation rate be much smaller than k B T '/%: the dynamics was therefore dubbed 
"quantum relaxational". Related ideas were applied to charge transport in the insightful 
article by Sondhi et al. J7|. Dynamic order parameter fluctuations also carry charge, and 
therefore inelastic collisions between thermally excited charge-carrying excitations will lead 
to a transport relaxation time r tr . As the typical energy exchanged in a collision is k B T, r tr 



is also of order t v , and therefore [[32 



1 k B T 

T^r h~' (L1) 

The missing coefficient in ( |1.1| ) is a universal number whose value will depend upon the 
precise definition of 1/V tr . It is perhaps worth noting explicitly here that (|1 . 1|) holds for 
all values of the dynamic exponent z. (A scenario under which ( |1 . 1| ) could be violated is 
discussed in Section |I V B|) . The collisions leading to ( |1.1|) conserve total energy and, in the 



case of translationally invariant systems, total momentum. Note that conservation of total 
momentum should not be confused with conservation of the total charge current — the latter 
is not conserved in any of the systems of interest here, even in the continuum scaling limit 
(the total charge density is, of course, always conserved). 

Now general considerations |33| suggest that there are two qualitatively different regimes 



of charge transport at non-zero frequencies 

• UT tr <C 1: the hydrodynamic, incoherent, collision-dominated regime, where charge 
transport is controlled by repeated, inelastic scatterings between pre-existing ther- 
mally, excited carriers; the conductivity should exhibit a 'Drude' peak as a function 
of frequency. 

• UTtr ^> 1: the high frequency, phase-coherent, collisionless regime, where the excita- 
tions created by the external perturbation are solely responsible for transport, and 
collisions with thermally excited carriers can be neglected. 

Essentially all previous theoretical analyses have been restricted to the collisionless regime. 

We can sharpen our description of the crossovers in the vicinity of uj ~ l/r tr by expressing 
them in universal scaling forms. We will consider systems where the order parameter quanta 
carry 'charge' Q. For the superfluid-insulator transition of bosons, we have Q = 2e and we 
measure charge transport by the dynamic conductivity a(u). For the Neel-paramagnet tran- 
sition in antiferromagnets, we have Q = g\i B (fi B is the Bohr magneton and g is gyromagnetic 
ratio) we and measure spin transport by the 'spin conductivity' (also denoted by u(uj)) which 
determines the spin current in response to a uniform gradient in an external magnetic field. 
In d spatial dimensions the dynamic conductivity a obeys at the quantum-critical coupling 

(1.2) 
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where z is the dynamic critical exponent, c is a non-universal microscopically determined 
quantity with the dimensions of (length) z (time) _1 (for z = 1, c is a velocity), and S(cJ) is a 
universal scaling function (note that we are using h rather than h to define the scale of the 
conductivity). We note again that the dependence on a universal function of huj/ksT, with 
no arbitrary frequency scale factors, holds for all values of the dynamic exponent z [|26 



Scaling as a function of huj/ksT was also noted recently by Sondhi et al. [Q], but they used 
it only to establish the frequency dependence of the conductivity away from the critical 
coupling, with an eye to understanding recent dynamical conductivity measurements in the 
quantum Hall system []33J (we will comment on these measurements in Section [IV D| ). Indeed, 
their scaling forms are consistent with a conductivity which is independent of huj/ksT right 
at the critical point, and hence with implicit or explicit assumptions in earlier theoretical 
results [|T4]-|25|. One of our new points here is that there is a non-trivial dependence on 
huj/ksT already at the critical coupling, and that this dependence means that previous 
analyses of the universal conductivity at the critical coupling either did not compute |H-F2T|, 
or were not particularly sensitive to |22|-f25|j, the value of the d.c. conductivity. 



Let us now turn to the expected dependence of the conductivity on the huj/ksT as 
expressed in the universal function £(u;). The discussion here will apply to realistic systems 
at their critical point for Tiuj and k^T much less than a non-critical microscopic energy 
scale e.g. the repulsion energy U between two bosons on the same site of a lattice model. 
The physical arguments below ( |1.1|) have been recast in a qualitative sketch of the real 
part of the function S(uJ) in Fig |l|. The d.c. conductivity is determined by the value of 
the real universal number S(0). At small uj there is a Drude-like peak coming from the 
energy-exchanging collisions among thermally excited carriers. At larger frequencies there 
is crossover to transport by particle-hole pairs created by the external source. As To — > oo 
we expect that S(cJ) ~ (— icu)( d ~ 2 ^ z so that a becomes independent of T in the collisionless 
regime. In d = 2, £(oo) is a real, finite, universal number determining the high frequency 
conductivity; S(oo) was the number computed in earlier analyses pi[-|2T|, and not the d.c. 



conductivity which is given by S(0). This is illustrated in Fig which plots the form of 
(h/Q 2 )a(uj,T — > 0) in d = 2: its value at uj = is given by S(0), while for all uj > it 
equals E(oo). It is likely, although not established, that E(0) > S(oo). Note the difference 
from Fermi liquid theory, where the Drude peak becomes a delta function with non-zero 
weight as T — > 0. In the present situation, the weight in the Drude-like peak vanishes like 
~ T( d + z - 2 )/ z as t — > 0, and in d = 2 and the Drude-like peak reduces to the single point 
uj = where the conductivity is given by S(0). The scaling properties of the d — 1, z — 1 
case are similar to those of a Fermi liquid. 

There is a certain critical phenomena/scaling perspective in which the result ( |1.2j ), and 
its implications discussed above, may seem quite natural, and even somewhat 'trivial' as they 
follow directly from the fact that uj and T have the same scaling dimension. Nevertheless, its 
importance has been overlooked in essentially all previous work. This is probably because 
there are complementary perspectives, more common among investigators in this field, from 
which ( |1.2j ) implies physics that is surprising and even somewhat radical. In particular, 
current ideas in quantum transport theory [^] and dissipative quantum mechanics H$6f can 



lead one to rather different picture, as we now itemize explicitly: 

• There have been a number of previous situations in which charge transport proper- 
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ties have been found to be universally related to the quantum unit of conductance, 
e 2 /h; these include the quantized Landauer conductance of ballistic transport in one- 
dimensional wires, and the universal conductance fluctuations of mesoscopic met- 
als |35| , |37|| . However in all previous cases, these universal properties have arisen in 
a "phase-coherent" regime, i.e. they are associated with physics at scales shorter than 
the mean distance between inelastic scattering events between the carriers. For the 
case of a d = 2 quantum critical point discussed above, the universal number £(oo) 
is associated with quantum coherent transport, and is therefore the analog of these 
earlier results. In contrast, the value of S(0) is controlled by repeated inelastic scat- 
tering events, and therefore the d.c. transport is clearly in what would traditionally be 
identified as the "incoherent" regime. Nevertheless, we have argued above that S(0) is 
a universal number, and remarkably, the d.c. conductance remains universally related 
to e 2 /h. 

The community has gained much intuition on the non-zero temperature transport 
properties of interacting quantum systems from recent studies of dissipative quantum 
mechanics and a number of related quantum- impurity problems |]36|j38| -^0| . The scaling 
properties of such models are given by the theory of boundary critical phenomema, 
in contrast to the bulk critical phenomena of interest in this paper. We discuss the 
transport properties of a class of boundary problems in Appendix the leading term 
in the transport coefficient is found to be independent of u/T, and dependence on u/T 
arises only upon consideration of subleading terms at low T. However, as we argue in 
Appendix |A], this behavior is understood by the fact that the fixed point controlling 
the low T behavior is simply a free field theory, and u/T dependence arises only upon 
considering the leading irrelevant operator. In contrast, the bulk theories of interest 
here differ in a crucial respect: they have interacting critical theories and therefore, 
we argue, contain u/T dependence already in the leading term, before the inclusion 
of any irrelevant operators. 

Scaling as a function of u/T does not hold for the Anderson localization transition of 
non-interacting electrons. The frequency obeys conventional scaling with dimension 
z, but the behavior of temperature is rather different and non- universal. Analytic 
theories for such transition are available only for d > 2, and in these the primary effect 
of temperature is in the non-universal, superlinear in T dependence of the phase- 
breaking rate which acts like a finite-size-like infrared cutoff to the critical properties. 



See also Section [TV B| . In contrast, the models of interest here have interactions, 
and inelastic scattering is central to understanding their universal critical properties; 
in some respects scaling with respect to temperature is simpler, as its naive scaling 
dimension of z is now valid. 

In this paper we will provide explicit results for the crossover function E(uJ) in a simple, 
disorder-free field theoretic model for the superfluid insulator transition that was introduced 
by Cha et al. ||15|| . Near the quantum-critical point, this model becomes equivalent to the 
familiar particle-hole symmetric 4 field theory with the effective imaginary time (r) action 

S = fJ kBT dr J d d x[\ [(^ Q ) 2 + c 2 (V x .0 Q ) 2 + K c + t o )^] + ^(0 2 ) 2 }. (1.3) 
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Here <p a is a n-component field and the action has 0(n) symmetry (the 0(n) index a is 
implicitly summer over). The spatial and temporal gradient terms are both second order, so 
the action has a "Lorentz" invariance with c the velocity of light, and as a result the dynamic 
critical exponent z — 1. The bare "mass" term has been written as m,Q C +to so that the T = 
quantum critical point is at to = 0, and uq measures the strength of the quartic non-linearity. 
The superfluid-insulator transition is described by the case n = 2 where ^ = <fti + i(f>2 is 
the usual complex superfluid order parameter. We shall also be interested in the case n = 3 



which applies to quantum-critical points in quantum antiferromagnets 2S 



It should be noted that the continuum model S for n — 2, with its double time derivative 



term, differs from the usual continuum action for non-relativistic bosons JO]] which has 
only a single time derivative of a complex scalar field. The total momentum and charge 
current of the latter model are proportional to each other. When these non-relativistic 
bosons are placed under the influence of an effective lattice potential, with an average of 
an integer number of bosons per lattice site, then after integrating out certain local high 



energy modes, one obtains the continuum model 5 as a low energy effective theory |42 



The charge current and momentum operators of S are now no longer simply related to 
each other as S has excitations with both positive and negative charges. As we will see in 
Section pTI] , the total momentum is proportional to the sum of the currents of the positive 
and negative charges, while the total charge current is proportional to their difference. This 
lack of a direct relationship between the charge current and the momentum should not seem 
surprising as lattice effects were required to obtain S, and are therefore implicitly accounted 
for in the continuum theory. 

It was asserted by Cha et al. that the T > transport properties of S were "patho- 
logical" in that there was zero resistance to charge transport in the continuum field theory in 
(|1.3[); they suggested that a non-zero resistance appeared only upon considering additional 
lattice corrections (beyond those required to derive the continuum theory S) , in which mo- 
mentum could jump in units of reciprocal lattice vectors at scattering events (the so-called 
"umklapp" scattering events). We shall show here that this is incorrect. The model S has 
a finite, and universal, d.c. resistance at any T > already in the scaling, continuum limit. 
Indeed, this is already clear from a recent study |43[] of S for the case d — 1 n = 3, where a 
simple and physically transparent argument obtained the exact (and finite) low temperature 
value of the spin diffusivity. Rather than being pathological, we claim that the transport 
properties of the disorder-free boson model S are generic, and essentially identical in their 
scaling structure to those of disordered boson systems. The error in Cha et al. |T5f appears 



to be due to their ignoring the presence of independent positive and negative charge excita- 
tions, and the resulting difference between the total momentum and the total charge current 
of S 0. 

As T appears in S only in the upper limit of the imaginary time integral, it is clear that 
the scaling form ( |1.2| ) is nothing but a standard finite-size scaling result for observables as 
a function of "wavevector" , u, and "inverse size", T. It might then seen that our job is 
relatively straightforward, and we merely have to obtain standard finite-size scaling results 
on the familiar model S. This is far from being the case. The point is that these standard 
results exist only at imaginary frequencies i2pnkBT/h (p > 1, integer), and we are especially 
interested in real frequencies -C ksT/h. More importantly, it has been shown p6| , p!3| , p9| , 
that the operations of analytic continuation and expansion in 1/n or e = 3 - d (which are 
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the only non-numeric tools for analyzing the critical point of S) do not commute. 

The proper tool for analyzing transport at the critical point of S is a quantum Boltzmann 
equation (QBE) for the charge carriers. In general, solution of such a QBE is a dauntingly 
difficult task, but we shall find that it is possible to reduce the solution to a simple linear 
integral equation in an expansion Jl9|,^9j in e = 3 — d (very similar results can be obtained in 
a related analysis of S in an expansion in 1/n [|15|,^(||13|] ; this will be described elsewhere, and 
we will only discuss the n = oo results here.). For small e, the non-linearities are weak, and 
it becomes possible to give a quasiparticle-like interpretation to the excitations of S at the 
quantum critical point. Such an interpretation will be useful in our intuitive understanding, 
and will help us use standard methods to simplify the QBE in the hu <C k^T limit. Although 
the quasiparticle interpretation fails for the physical case e = 1, we do not expect any 
qualitative change in the structure of our results for larger e. The QBE formulation is 
quite general, and the quasiparticle representation is mainly a useful technical tool towards 
obtaining its numerical solution. In a more general context, our approach may considered as 
an expansion in powers of the anomalous dimension, 77, which is responsible for replacing the 
quasi-particle pole by a continuum at the T = critical point. The continuum associated 
with a non-zero 77 is important only for tkv ^> ksT, while for fno ~ ksT very different 
thermal damping processes quench the critical fluctuations [ 2S| , and are best treated by the 
QBE. This thermal damping also acts as an effective infrared cutoff which ensures that no 
qualitatively new physics emerges at higher orders in the expansion in powers of r\. 

The basic structure of the dynamical conductivity at the critical point of S for small e 
is illustrated in Fig || Notice the presence of "boundary layers" which make the analysis 
of the e — > limit quite subtle, and is responsible for the non-commutativity of analytic 
continuation and the naive e expansion noted earlier. The hydrodynamic regime of the 
conductivity (denoted later in the paper by 07) is represented by a Drude peak of width in 
frequency uj ~ e 2 ksT '/h and has a height of order 1/e 2 . In particular, the d.c. conductivity 
is determined by the universal number S(0), for which we find for n = 2 

m = ™» (1.4) 

to leading order in e; the structure of the higher-order corrections to ( |1.4p is quite complex, 
and was generally discussed in Ref [ 29| . Determination of the coefficient in (|1 .4\ ) required the 



numerical solution of a QBE, and the uncertainty in the numerics is believed to be restricted 
to the fourth decimal place. In d = 2, Q1.4Q implies a universal d.c. conductivity 

4e 2 4e 2 

o-(O) = 2ttS(0)— « 1.037— (1.5) 
h h 

at the superfluid-insulator transition. This result is remarkably close to the self-dual value 
4e 2 /h PByTTI , and to the results of many experiments 0; we will comment further on this 
in Section [IV E . 

There is a clean separation between the hydrodynamic and collisionless regimes of o: 
the latter does not begin until uj ~ c 1//2 /cbT/^ (Fig |3|). The T = collisionless transport is 
characterized by S(cJ — ► 00) for which []19| we have at n = 2, 3 
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As noted earlier, for e = 1 (d = 2) E(oo) is a pure number, but notice that it bears no 
relationship to X(0); indeed E(oo) is of order unity as e — > 0, and so these two quantities 
are of distinct orders in e. For the superfluid-insulator transition in d = 2, fll.6| ) gives a high 
frequency conductivity of 0.3927 x (4e 2 /h), which is in rough agreement with other analyses 



in the collisionless regime (19 . 

The body of the remainder of the paper is devoted to obtaining the above properties 
of the model S. Readers not interested in calculations specific to the model S should now 
go directly to the concluding Section [IV] where we discuss implications of our results for 
a number of experimental systems. In Section [TT] we will obtain one loop results for the 
transport properties of S using the familiar Kubo formalism. Then Section |TTT| will include 
two-loop effects using a quantum transport analysis needed to describe the hydrodynamic 
regime, and obtain ( |1.4| ). 

In all of Sections |I|, [ITJ, an d the appendices we will work in units in which h = ks = c = 1; 
we will reinstate these constants in Section |IV]. 



II. ONE LOOP RESULTS FROM THE KUBO FORMULA 



We will begin our analysis of the transport properties of S by examining the results 



of a direct evaluation from the Kubo formula [|X~5] , [X~9|] , but now working at T > 0. The 
physical interpretation of the results will motivate an analysis using a quantum Boltzmann 
equation p3| , ^6| , ^7| which will be carried out in subsequent sections. 

The standard Kubo formula relates the conductivity to a two-point correlator of the 
conserved 0(n) current. We introduce an external vector potential A associated with the 
0(n) generator which rotates <ft a in the 1,2 plane; the spatial gradient term in S then 
undergoes the mapping 



(V,0i - QA0 2 ) 2 + (V,0 2 + QA0O 2 + £ (V a 



(2.1) 



a=l 



«=3 



The associated 0{n) current is then 5S/5A. We evaluate its two-point correlator using 
standard diagrammatic perturbation theory to first order in uq in the expansion in e, or to 
leading order in the large n expansion (in which uq ~ The first order vertex correction 

vanishes because the interaction is momentum independent and the result in both cases is 
given simply by p~5|JT9[] (recall we are using units here in which U — ks — c — 1) 



2Q' 



tJ2 



d d k 

(27T) a 



2kl 



(e^ + k 2 + m 2 )((e n + uj n ) 2 + k 2 + m 2 ) + k 2 + m 2 



(2.2) 



The first term is the 'paramagnetic' contribution, while the second is the 'diamagnetic' 
term proportional to the density |48]. Here e n , cu n are Matsubara frequencies, k x is the 
x component of the d dimensional momentum k, and k = |k|. The "mass" m in the 
propagators is computed in Appendix (where it is referred to as m(T)) using the e = 
3 — d expansion developed in Ref p9| , and depends universally upon T and an energy scale 



8 



measuring the deviation of the ground state from the critical ground state. At the critical 
coupling, to leading order in e, 



m 



n + 2\ 2vr 2 T 2 



n 



at tr 







(2.3) 



The result Q2.2|) holds for all n, and there is no n dependence at this order in the e expansion, 
other than that through m. The large n expansion has an identical structure at n = oo, the 
only difference being in the value of m. Detailed universal expressions for m were given in 
Ref. |l3j in d — 2; at the critical coupling, the analog of ( |2.3|) at n = oo is 



m 



2 In 



'V5 + l' 



T 



at tr 



0, d 



(2.4) 



The remaining analysis of this section will apply both to the e and 1/n expansions, the only 
difference being in the values of m given above. 

Now insert 1 = dk x /dk x in front of the diamagnetic term in (|2.2j ) and integrate by parts. 
The surface terms vanish in dimensional or lattice regular izat ion, and the expression for the 
conductivity becomes 



a[lLOr, 



2Q' 



OJr. 



tJ2 



d d k 



2kl 



(2Ti) d e 2 n + k 2 + m 2 



1 



1 



(e n + uj n ) 2 + k 2 + m 2 e 2 +k 2 + m 2 



■ (2-5) 



We now evaluate the summation over Matsubara frequencies, analytically continue to real 
frequencies. The resulting <j(uj) is complex, and we decompose it into its real and imaginary 
parts cr{ijj) = cr'(^) + ia"{uj). We will only present results for the real part ct'(uj), and the 
imaginary part <?"{uj) can be obtained via the standard dispersion relation. 

We find that the result for a'iuj) has two distinct contributions of very different physical 
origin. We separate these by writing 



cr'(uj) = c7j(lu) + o ! h (uj). 
The first part, (t'j(u), is a delta function at zero frequency: 

d d k k 2 x ( dn{e k ) 



a'j{uo) = 2ttQ 2 S(u; 
where n(e) is the Bose function 

n(s) 



(2k)* 4 



1 



e e/T _ 1 ' 

and the excitations have the energy momentum relation 

e k = \/k 2 + m 2 



(2.6) 



(2.7) 



(2.8) 



(2.9) 



We will discuss the physical meaning of the delta function in ( f2.7Q below, and obtain a 
separate and more physical derivation of its weight in Section [III A| . The second part, 
cT' n (uj) is a continuum above a threshold frequency of 2m: 
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(/„(*,) = vrQ 2 / 7^^(1 + MskMH - 2e k ) 

= ^pfyM - 2m) ^ 2 ~J m y /2 [1 + 2n( W /2)] M d " 2 , (2.10) 

where Sd = 2/(r(<i/2)(47r) d//2 ) is a standard phase space factor. 

At the critical point, to = 0, it can be verified that the above results for cr(u) obey the 
scaling form (|1.2j ) with z = 1; explicit results for the function E(u7) will appear below and 
are sketched in Fig |j. 

We now discuss the physical and scaling properties of the two components of the con- 
ductivity in turn: 



1. aj 



This is a zero frequency delta function, and is present only for T > 0. Clearly, it must be 
interpreted as the contribution of thermal excitations which propagate ballistically. Indeed, 
to first order in e [29|, or at n = oo [13|, the excitations are simply undamped particles 
and holes with an infinite lifetime and energy momentum relation It is necessary to 
go to second order in e, or to first order in 1/n, to include collisions which will give the 
quasiparticles a finite lifetime. We will show in Section [TTJ that these collisions also broaden 
the delta function in 07. The magnitude of the broadening is expected to be determined by 
the inverse lifetime of the quasiparticles; at the critical point, this inverse lifetime is of order 
e 2 T [^] in the e expansion, or of order T/n [13] in the large n theory. The typical energy of 
a quasiparticle at the critical point is of order T, and so the quasiparticles are well-defined, 
at least within the e or 1/n expansion. Notice, however, that the quasiparticle interpretation 
breaks down at the physically important value of e = 1, n = 2. The discussion in Section [Tl I 
will take place within the context of the e expansion, and will use the quasiparticle interpre- 



tation and the arsenal of powerful techniques available |33| , ^6| , ^7| to describe their quantum 
transport. 

The expression ( |2.7| ) is valid everywhere in the normal phase, but here we evaluate it 
explicitly only at to = 0. Consider first the e expansion. The coefficient of the delta function 
is a function of the ratio m/T, but notice from ( |2.3j ) that m -C T at for small e. Evaluating 
in this limit we find for e small 



a'j(oj) = 2nQ 2 T d - 1 5(uj) 
= 2irQ 2 T d - 1 5(uj) 



m 



18 8ttT 



+ • 



1 ^/2(n + 2) 



1/2 



8 \3(n + 8) 



+ 



(2.11) 



Actually the expression ( |2.7| ) is good to order e but we have refrained from displaying the 
next term as it is rather lengthy. The first term in (|2.11| ) is obtained by evaluating (|2.7p at 
m = 0, d = 3; the second term is from an integral dominated by small ~ m « T and 
hence the Bose function can be replaced by its classical limit. An important point to note is 
that the current carried by the thermally excited carriers is dominated in the leading term 
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of (|2.11|) by momenta k ~ T ^> m. This will be useful to us in the analysis of collisions in 
Section |T| where we will simply be able to set m = to obtain the leading term. 

Turning to the large n theory, we evaluate ( |2.7| ) in d = 2 using the value of m in ( p.4| ). 
and obtain after an integration by parts and rescaling of variables 



Q 2 T 
2 

Q 2 T 




6(u) 



S(u) x 0.689402548116632 . . . 



(2.12) 



where 



ory 



13,49 



21n((v5 + l)/2) is a number which plays a central role in the large n the- 
Notice that as m ~ T, we have now been unable to approximate Ek ~ k to get 
the leading result, as was done in the e expansion. 

An interesting numerical property of the above results in d = 2 is worth noting explic- 
itly. The spectral weight of the delta function to leading order in the e expansion is, from 
( [2. 11| ), Q 2 T times tt/9 = 0.3491 . . . (recall that this number was obtained by evaluating a 
momentum space integral in d — 3). The same quantity at n — oo, from (|2.12|) , is Q 2 T 
times 0.3447 . . . (obtained now by evaluating a different momentum space integral in d — 2), 
which is remarkably close. Later in Section we will consider broadening of the delta 
function in the e expansion, and we will work in the approximation in which the spectral 
weight is Q 2 T times 7r/9 (see Eqn (|3.34j ) later); the present numerical "coincidence" suggests 
that the numerical values of the leading order e result are quite accurate. 



2. a n 

This is the continuum contribution to a which vanishes for uj < 2m. At this order in 
e (1/^) there is a sharp threshold at uj = 2m but we expect that this singularity will be 
rounded out when collisions are included at order e 2 (1/^): we will not describe this rounding 
out it in this paper, however. Although they have a strong effect at the threshold, collisions 
are not expected to significantly modify the form of u' n {uj) at higher frequencies where the 
transport is predominantly collisionless. In particular, the uj — ► oo limit is precisely the 
T = result obtained earlier 



a n (u -> oo) = -yq-W (2.13) 



III. QUANTUM TRANSPORT EQUATIONS 

The general analysis of higher order corrections to a is quite complex, and so we will 
confine ourselves in this section to the answer to a single question: how does the S(u) term 
in a'jiui) broaden ? 

We will address this question exclusively in the context of the e expansion, and gen- 
eralization to the 1/n expansion will be discussed elsewhere. In principle, the answer can 
be obtained by including the 0(e 2 ) correction to the self energy and accounting for the 
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associated infinite-ladder vertex corrections |48| . However this method is quite inconvenient 
as it does not allow for easy separation of the distinct phenomena in different frequency 
regimes. Instead, we shall use the (in principle) equivalent [E8| quantum transport formal- 
ism of Kadanoff and Baym |33|j46|j47|| . The physical content of this formalism is transparent 
at all stages, and the approximations necessary to focus on the low frequency conductivity 
are readily apparent. In particular, we can drop the terms leading to <Jh{oj) at an early 
stage. 

The transport equation is best studied in the Hamiltonian formalism by casting it in 
terms of the weakly interacting "normal modes" . So, we begin by writing down the Hamil- 
tonian associated with S (reminder — we are using units in which h = ks = c = 1): 



7i = 7~Cq + Ti. 



int 



n 



ext. 



The first term, Tio is the free particle part of S 

1 



Ho 



*l + (V x (j) a ) 2 + m 2 (j) 2 



(3.1) 



(3.2) 



where 7r a (x, t) (t is real time) is the canonically conjugate momentum to the quantum field 
4> a (x,t) and satisfies the equal-time commutation relations 



(x,t),7T /3 (x',i)] =z5 a/3 5 d (x-x'). 



(3.3) 



For future convenience, we have already included the Hartree-Fock correction from the 
interactions into the 'mass' m 2 (these correspond to self-consistently summing all the one- 
loop tadpole diagrams, as discussed in Appendix |D|). The second term, TCi nt is the quartic 
interaction 



n 



int 



Uo 

4! 



d d x( 



l2 \2 



(3.4) 



and it is understood that the Hartree-Fock term arising from 7i int will be omitted. Finally, 
7i ex t contains the coupling to the external space and time- dependent potentials f/ a (x, t) 
(o= l...n(n-l)/2) 



H ext = Qj d d xU a (x,t)L a a ^ a (x,t)M^t). 



(3.5) 



Here the L a are n x n real, antisymmetric matrices that are i times the generators of the 
Lie algebra of 0{n). The U a are coupled to the conserved 0{n) charge densities ]50|] of TC 
+ TYint- We shall be interested only in the linear response of the current to the 'electric 
field' E a = — V x C/ a (x, t), and it will be assumed below that E a is independent of x. Notice 
that we are making a gauge choice different from that in ( |2.1|) , and coupling now to the 
scalar and not the vector potential; this is for convenience, and should not change the final 



gauge-invariant results. The 'charge' current J a is defined by the expectation value [50 



(3.6) 



and it will also be independent of x. Making the Fourier expansion 
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E a (t) = / ^E»e 



-iuit 



(3.7) 



and similarly for J a , we can define the dynamical conductivity, cr(uS) by the expected linear 
response relation 



J a M = <xME a M 



(3.8) 



For completeness, let us also note here the expression for the total momentum [50], P 

P = (n a V x 4> a ) . (3.9) 

Notice that it is quite distinct from 3 a . In particular, in the absence of an external potential, 
P is conserved (i.e. it obeys an equation of the form dtP + V- T= for some T), while 3 a 
is not. 

We now make the mode expansion 



d d k 

(2n) d ^2i k 
d d k m 



0a(x,t) = / ^-^-^=(a Q (k,t)e ikx + 4(k,t)e 

7T Q (x,t) = -% 



ikx 



(2n 



(/V ^(a a (k,t)e tkx -at(k,t)e- k - 



(3.10) 



where the a(k, t) operators satisfy the equal-time commutation relations 

a Q (k,t),4(k',t)] =6 aP (2n) d 6 d (k-k') 
[a a (k,t),ap(k',t)] =0. 

It can now be verified that ( |3.3| ) is satisfied, and Ho is given by 

d d k 



7~tc, 



(27T) 



4(k,t)a Q (k, t) + 1/2 



(3.11) 



(3.12) 



We will also need the expression for the current J a in terms of the a and a) . We will only 
be interested in the case where the system carries a position-independent current: for this 
case, inserting ( |3.10| ) into ( |3.6| ), we find 

r(t) = jj(f) + j?,(t) 

c/ d fc k 



__(ot(k,t) a/3 (k,t) 
r u {t) = -iQL a ap j 



{2ir) d e k 
d d k k 



{27r) d 2e k 



-k,t)al{k,t)\ + H.c. 



(3.13) 



It should be evident that processes contributing to J" 7 require a minimum frequency of 2m, 
and so i a H only contributes to an(u). We will therefore drop the Jjj contribution below 
and approximate 3 a ~ J". The ease with which the high frequency components of &(oj) 
can be separated out is an important advantage of the present formulation of the quantum 
transport equations. 
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A central object in transport theory is the Green's function |33| , [47 
d d K 



g<Jk,n,R,t) 



(2tt) c 



aUk - K/2, t - ti/2)a a (k + K/2, t + h/2) 



(3.14) 



For the case of a system carrying a spatially independent current, g < will be independent of 
R and this will implicitly be assumed below by dropping the R argument. We also define 
the particle distribution function 



f a p(k,t) = J — g <p(k,n,t) 

in terms of which the current is 

r{t) = iQL a aP J 



d d k k 



/a/?(M) 



(2ir) d e k - 

The corresponding expression for the momentum density is 

d d k 



P(f) 



]kf aa (k, t). 



(3.15) 



(3.16) 



(3.17) 



Notice the difference in the structure of the 0(n) indices between (|3.16|) and ( |3.17 ). 

We note here that our formulation of the transport theory in terms of Green's functions 
of the a a , a^ a rather than those of the (j) a , n a was motivated in part by recent exact results 
in d — 1 for time- dependent correlations of models equivalent to S for n — 1,3 |]5T| , ^3| . In 
the latter cases it was evident that the physics is most simply described by following the 
propagators of particles created/annihilated by aj,,/a a through their collisions. 



A. Collisionless transport 

In this section we will examine the transport equations for the collisionless case where 
'Hint — 0; we remind the reader that interactions have already been included at the Hartree- 
Fock level in 7io (see Appendix 0). Strictly speaking, we also have to remember that the 
mass m can in general depend upon E a ; however for the case of a momentum- independent 
local interaction uq, such a 'vertex' correction vanishes to lowest order in E a , and will 
therefore be omitted from our discussion. 

While it is possible to discuss the general 0{n) case, in the interest of simplicity and to 
keep the physical content transparent, we will restrict our attention here to the special case 
n = 2. The generalization to n > 2 is discussed in Appendix ||. 

For n — 2 (a = 1,2), there is only one real antisymmetric matrix, and therefore the 
index a can be dropped. We choose L 12 = — L 2 ,i = 1 and L^x = L 2< 2 = 0. This matrix is 
off-diagonal and it is helpful to transform to a basis where the external field is diagonal. We 
therefore define 

a ± (k,t)^ ai(kt)± ; a2(M) (3.18) 
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The current now becomes 

d d k ^ . k 



J 



Q J 7^E A -(4(k,t)a A (k,t) 



-Qj^^UKt) (3.19) 

where the index A is assumed here and below to extend over the values ±1, and fx are the 
particle distribution functions which are now diagonal in A space. Notice that there are 
two species of charged particles with charges ±Q: these are the particle-like and hole-like 
excitations of the bosonic insulator. Let us also note the expression for the momentum 
density 

/rjdh. 
(^X>WM) (3.20) 



An important difference between ( |3.19| ) and ( |3.20|) is the A inside the summation in ( 3.19|) 



which is absent from ( ft.20| ). Thus the 'charge' current is proportional to the difference of 
the particle and hole number currents, while the momentum density is proportional to their 
sum. 



It is now easy to use standard methods to derive the following transport equation 

in the collisionless limit described earlier 

(1 + X QV(t) ■ J-j / A (k, t) = 0. (3.21) 

In deriving this equation we have made approximations to the charge density appearing in 
Hext similar to those made for J: upon expressing 7i ex t in terms of the a, a' we have dropped 
all terms involving the product of two a's or a''s as these will only contribute to the high 
frequency an- The equations (|3.19| , [3.2il ) are therefore accurate to first order in uq provided 
we are limiting ourselves to frequencies u <C 2m. 

Let us now solve (|3.19|J3.2T| ) in linear response. In the absence of E, the distribution 
function has the equilibrium value given by Bose function / A (k, t) = n{ek). We Fourier 
transform from time, t, to frequency u, and parametrize to linear order in E: 

/ A (k, uj) = 2n5(uj)n(£ k ) + AQk ■ E(u)ip(k, to), (3.22) 

where we have used the fact that only E breaks spatial rotation invariance and 0(2) sym- 
metry to conclude that ip is independent of lz/k and A. Now inserting in ( fj.21|) , and using 
dek/dh = h/ek it is simple to solve for ip to leading order in E: 

= _|_L(-«g£!)) (3.23) 
-lUJE k \ ds k J 



Finally we insert in ( |3.19| ) and deduce the conductivity 



2Q 2 f d d k hi ( dn(e k )\ 
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The real part of this agrees with (|2.7|) . Notice that the leading factor of 2 comes from the 
sum over A. The current is therefore carried equally by the thermally excited particles and 
holes: they move in opposite directions to create a state with vanishing momentum but 
non-zero charge current. 

We will see in the next section that this charge current can be relaxed by collisions among 
the particles and holes. There is no need to invoke umklapp scattering as a momentum 
sink/source (as was done by Cha et al. |T5|) as the current carrying state does not have a 
net momentum to begin with. 



B. Collision-dominated transport 

We have seen in Section [111 A| that to order u (e) the transport is described by the 
ballistic motion of undamped particles of two charges. We now consider the collisions of 
these particles which appear at order e 2 . As we noted in the discussion below fl2.11| ) the 



typical energy of a particle contributing to the transport was of order T. Their collisions 
will lead to a broadening of the single quasiparticle pole of order e 2 T [^UJ; for e small, 



this broadening is weak, and it is then permissible to argue in terms of a quasiparticle 
interpretation. 

Applying the standard methods of transport theory discussed in Refs [B3,l|7[], we gen- 



eralize (|3.21| ) to include the collision term which appears at order Uq, and obtain to linear 
order in E: 

(d_ d_\ . , _2u| r d d h d d k 2 d d k 3 1 

\dt + ^ ' dk) M ' j " 9 J (2ir) d (2n) d (2n) d I6s k e kl e k2 e ks 

x (27r) d 5(k + k a - k 2 - k 3 )2vr5(e fc + e kl - e k2 - e ks ) { 
2/ A (k, t)/- A (ki, t) [1 + / A (k 2 , t)} [1 + /_ A (k 3 , t)} 
+ / A (k, t)f x (k u t) [1 + / A (k 2 , t)} [1 + / A (k 3) t)} 

- 2[1 + / A (k, t)\ [1 + /_ A (k 1; t)]/ A (k 2 , *)/_ A (k3, t) 

- [1 + / A (k,t)][l + / A (k 1 ,t)]/ A (k 2 ,t)/ A (k 3 ,t)}. (3.25) 



The collision term on the right-hand-side of (|3.25|) can also be obtained by a simple argument 



based on Fermi's golden rule of the type described in Ref. pB|. A number of simplifications 



have been made in deriving ( |3.25| ), and we now describe and justify them: 

• The collision term is initially expressed in terms of full two-point Green's functions like 

However, as is conventional p3|JT7[1 , we assume we can neglect damping in these 
Green's functions and express them in terms of the particle distribution functions f a p. 
This is permissible for a system with well-defined quasiparticles, as is the case here for 
small e. This approximation will not affect the conductivity at order e 2 . 

• In addition to the collision terms, there are also self-energy terms affecting the quasi- 
particle energies at order Uq. These include terms that couple the usual Green's func- 
tions to the anomalous ones (involving expectation values of pairs of creation or anni- 
hilation operators). Such terms will, in general, modify the conductivity at order e 2 . 
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However, they do not affect the nature of broadening of the delta function in <7j(u) 
The total spectral weight in cr'j(u) will change from that in ( 2.11| ) at order e 2 due to 



these terms, but the functional form of of <j'j(uj) at frequencies of order e 2 T will not be 
affected. 

We have neglected collisions which involve creation or annihilation of particle-hole 
pairs as they have negligible phase space. Thus a collision in which e.g. a positively 
charged particle of momentum k turns into two positively charged particles and a 
negatively charged hole with momenta ki, k2, and k3 respectively is permitted by 
the symmetries of the problem. However, it remains to evaluate the phase space over 
which such collisions conserve total energy and momentum. Notice that the 'mass' m 
of the particles/holes is of order yfeT (Eqn (|2.3|) ) while their momentum is of order T. 
So to leading order in e we may just replace the energy momentum relation ( |2.9| ) by 
(see also the discussion below (|2.11|) ) 



e k = k. (3.26) 

We will use this simplified energy-momentum relation throughout this subsection. 
The particle-hole pair-creation collision now requires that k = kx + k 2 + k 3 and k = 
k\ + k 2 + k 3 . This is only possible if all three momenta are collinear, and this process 
therefore has vanishing phase space. More generally, for a non-zero m, the phase space 
vanishes as e — > 0. 



We now insert the parametrization (ft.22 ) in ( 3.25Q , linearize the resulting equation in 



the external electric field E, and obtain a linear integral equation for the unknown function 
ip(k, u). Further, to leading order in e, we may set d = 3 in the collision term in (|3.25|) , and 
replace the coupling w by its fixed point value for n = 2 (see Appendix |B|) 

247r 2 e 

uo = — — . 3.27 

5 

Further details of this linearization procedure are given in Appendix y. The final integral 
equation for ip(k,u) can be written as 



1 dn(k) , 



oo 



-icoip(k,uj) + -—^ = -e 2 I dk 1 [F 1 (k,k 1 )ip(k,uj)+F 2 (k,k 1 )ip(k 1 ,co)). (3.28) 



o 



The expressions for the functions F 12 are quite lengthy and are discussed in Appendix y. 

It is now useful to scale out the dependence of all functions on e and T so that the final 
integral equation is written in dimensionless form, and has all couplings of order unity. From 
the expressions in Appendix we know that -F 12 are homogeneous functions of momenta 
and T which satisfy 

F i (k,k 1 ) = $ i (k/T,h 1 /T) (3.29) 

for some $j, with i = 1,2. (Here, and everywhere in this section, the dimensional analysis 
is performed for at d — 3: it is not difficult to extend it to general d, but we will not in the 
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interests of simplicity, and because it is unnecessary to obtain results to leading order in e.) 
Now by examining ( p. 28 ) we see that it is useful to introduce the function \1/ defined by 



if>(k,w) 



2^3 



T' e 2 T 



(3.30) 



In terms of $j, $ the integral equation fl3.28j ) takes the dimensionless form 

1 



— iu^f(k, u>) — 



Aksmh z (k/2) 



dfci k 1 )'^(k, uj) + <& 2 (k, k^i^kx, uj) , (3.31) 



where k = k/T and to = uj/e 2 
a i [yS] obeys 



uj/e T. Now, from ( |3.19|) we see that the conductivity 



2rpd-2 



Q T 



where the scaling function Ej is given by 



UJ 



e 2 T 



3tt 2 



'—3 — — 

k dk^(k,u)) 



(3.32) 



(3.33) 



We can now already see from the structure of (|3.32|) that <j'j(uj = 0) has a value of order 
T d ~ 2 /e 2 , and that the width of the 'Drude' peak will be of order to ~ e 2 T. Further, the 
collision term will not modify the total spectral weight in 07, which will therefore be identical 
to that in ( |2.11| ); this implies that the function £/, as defined by ( |3.31p and ( |3.33| ) should 
satisfy 



Jo d ^^ = T8- 



(3.34) 



It should be noted that this sum rule is special to the leading order in e being considered 
here. For e of order unity, there is no sharp distinction between 07 and ajj and there is no 
sum rule: indeed the integral in ( |3.34| ) when carried out over the total a will be divergent. 
For any realistic lattice model there is a large microscopic energy scale (~ U the repulsion 
between bosons on the same site) beyond which the universal scaling results do not apply, 
and the entire spectral weight (including frequencies beyond U) is not divergent; this latter 
spectral weight satisfies a sum rule related to non-universal microscopic quantities, and is 
unrelated to the universal result ( |3.34j ). 

It now remains to numerically solve (|3.31|) to determine \l/(/c,cD), and then to obtain 
from ( |3.33| ). The integral equation ( |3.31| ) was solved by a straightforward numerical 
iteration, and we found very rapid convergence to its unique solution. We show a plot of its 
solution at a few values of u in Fig The final result for the universal function (uj) is 
shown in Fig |[ 



IV. CONCLUSIONS 

The central message of this paper is simply stated: understanding the universal d.c. 
conductivity of a two-dimensional system at its quantum critical point requires a non-zero 
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temperature analysis of the hydro dynamic, collision-dominated regime where %lo <C k B T. 
The transport in this low frequency regime is incoherent, but nevertheless the remarkable fact 
is that the conductivity is still a universal number times e 2 /h. Computations carried out 
exactly at T = 0, with u — > 0, do not yield the d.c. conductivity, and are controlled by very 
different physical processes involving phase-coherent, collisionless transport. A related com- 
ment is that a theoretical analysis for the d.c. conductivity must necessarily be formulated 
in real time, as the imaginary time Matsubara frequencies are of order 27r/csT/fi or larger, 
and cannot easily capture the singular structure in Fig |2|. These criticisms apply not only 
to computations of the universal conductivity at the superfluid-insulator transition JTjj-|25[] 
but also to the transitions between the Hall plateaus |5^J53| , all of whom computed the 
analog of S(oo). There are also mappings between these models |54|, but they presumably 
hold for all values of the ratio huj/ksT. There are also special self-dual models |55| , |56| in 
which the conductivity is claimed |57]] to be independent of hu/k B T, but this is not expected 
to be the generic physical situation. The computation of the finite T conductivity at the 
superfluid-insulator transition in one dimension by Giamarchi and Schulz |58|| (and close to 
one dimension by Herbut ]59|]) seems closer in spirit to our approach, but it would be useful 
to have an explicit computation as a function of hu/ksT to verify this. 

Many experiments on two-dimensional films appear not to observe a universal conduc- 
tivity at the superfluid-insulator transition. Our results imply that this is very likely due 
to crossovers caused by inelastic scattering mechanisms other than those contained in the 
critical theory. Measurements of the conductivity as a function of frequency (Section [IV D| ) 
should help disentangle these effects. 

In the remainder of this section we comment on some general experimental and theoretical 
issues related to the results of the paper. 



A. Imaginary time Monte Carlo 

Most existing numerical studies [1^,22-25] of a at a two-dimensional quantum critical 
point used a Monte Carlo simulation in imaginary time. The simulation measures the 
values of o at the non-zero Matsubara frequencies i.e. at a(2imTi) where n is a non- 
zero, positive integer. The limit Tl(z, \z\ — > oo) is expected to be the same number S(oo) 
for all values of the phase of z, and therefore its value can be deduced from a relatively 
straightforward analysis of the numerical data. The value of S(0) is more problematical. 
All of the interesting structure in X(z) discussed in this paper occurs for z of order, or 
less than unity, and it is difficult to extract this information from its values at the nonzero 
quantized Matsubara frequencies. Cha et al. numerically examined the model S and 



found little dependence on LO n /T. Note that this is the model for which we have computed 
the strongly uj/T dependent conductivity here, but, as expected, this is apparently not 
evident at the imaginary Matsubara frequencies. Similarly, Wallin et al. ||25|| , also found 
little evidence for a significant dependence of the critical conductivity of disordered models 
on u n /T. 

We are of the opinion that it will be difficult to determine S(0) from this method, unless 
highly accurate numerical results are obtained at the imaginary frequencies. The difficulty 
is also apparent by a glance at Fig |2|: at T = 0, the value of S(0) appears only at a single 
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point which carries zero weight under any integral over frequencies. The accurate numerical 
data should then be analyzed in the following manner. First, from observations at a number 
of different values of u n and T, the universal scaling part of a, dependent only on the ratio 
u) n /T, should be obtained. Note that this universal scaling part reaches a constant value 
((4e 2 //z)E(oo)) as w n — ► oo, while the full a vanishes as for frequencies larger than 

a microscopic lattice scale. Then a Pade analysis, general enough to allow the structure 
in Fig [I], should be used to analytically continue only the universal scaling part to real 
frequencies. 

It would also be interesting to explore newer methods: perhaps examining an open 
system in which it is possible to have a net current in thermodynamic equilibrium, or doing 
a computation in real time. The exact diagonalization approach of Runge ^Tf can perhaps 
be extended to T > without a great deal of difficulty 



B. Dangerously irrelevant interactions 



It is possible to violate the basic scaling result ( |1 . 1| ) if quantum-mechanical interactions 



between the elementary critical excitations happen to be dangerously irrelevant |K| ; in that 
case we expect 



u rpl+0u/Z 



Ttr 



(4.1) 



where u is proportional to some power of the irrelevant interaction, 9 U > is the asso- 
ciated crossover exponent, and z is the dynamic critical exponent. Anderson localization 
transitions, with interactions leading only to a infrared cutoff as a phase-breaking rate, are 
a realization of such a scenario. (Quantum-impurity critical points |36|j38| -|4"0"| have some 
similarities to bulk systems with dangerously irrelevant interactions, and are discussed in 
more detail in Appendix However, this scenario is much less likely to be realized at bulk 
two-dimensional quantum critical points, and we consider it unlikely that interactions can be 
neglected for the superfluid-insulator transition. The quantum Hall transition has primarily 
been studied using non-interacting electrons [61] , but there is evidence that Coulomb inter- 



actions are relevant at quantum Hall transitions [62.63] . Experimental measurements of z 
in the quantum Hall system [64]] indicate the value z = 1 for the dynamic critical exponent, 
and this is incompatible with free electron models. The a.c. conductivity measurements of 
Engel et al. [[34] in a quantum Hall system (discussed further in Section [IV D| below) show a 
characteristic frequency scale u ~ ksT/h which is inconsistent with the result QL1| ) for the 
case of irrelevant interactions, but is consistent with the interacting theory result 



C. Luttinger Liquids 

A well studied critical system in d = 1 is the gapless Luttinger liquid ground state of 
interacting fermions or bosons away from commensurate filling fractions. Here we discuss 
how this familiar system fits into the general framework of this paper. The non-zero T 
conductivity of the critical theory of the Luttinger liquid is given by 
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(4.2) 



where K is some T independent constant. Notice that there is no broadening of the delta 
function, even for T > 0, in the scaling limit. (Related considerations also apply to Fermi 



liquids in d > 1 |65| , |66| , but for many purposes these are better thought of as the analog of 



a\u) = K(k B T)~ 1 5 — - . (4.3) 



Goldstone phases rather than critical phases.) Let us now rewrite (|4.2|) as 

Mr, 

Now notice that (|4.3|) is consistent with the general scaling form ( |1.2| ) for d — 1, and the 
known dynamic exponent z — 1. The scaling function £ takes the simple and singular form 
S'(oJ) = 5(uJ). For this special form of S, the limits u —>■ and T — > do happen to 
commute. 

These singular properties of the critical theory of the Luttinger liquid are clearly a con- 
sequence of the absence of scattering between carriers in the critical theory: like many other 
critical theories in d — 1, it is conformally invariant, and correlation functions (including 
those for T > 0) factorize into independent components given by the left and right movers. 
In contrast, for critical theories in d > 1, like the one studied in this paper, no such analogous 
factorization exists. Further, for d > 1 models below their upper critical dimension, there is 
scattering between the carriers already in the scaling limit. To introduce scattering in the 
Luttinger liquid, it is necessary to go beyond the scaling limit, and consider corrections to 
scaling |[45|| . There corrections are therefore dangerously irrelevant, and will destroy simple 
delta function form of the conductivity in Q4.2|) , and the limits to — > and T — > will no 
longer commute |[45||. 



D. Measurements of the a.c. conductivity 

A finite u measurement of the conductivity in a situation related to that discussed in 
this paper has been carried out by Engel et al. [j34|| . They examined the transition between 
integer quantum Hall plateaus by studying the dependence of the conductivity on lu, and 
the deviation of the field from its critical value 5 = (B — B c )/B c . This should obey |7j the 
scaling form 

e 2 ~ / huj 5 \ 



which generalizes ( |1.2| ) to 5 ^ 0. In their analysis, Engel et al. J34j focussed mainly on the 8 
dependence at Tioj ^> A^T: they measured the width of the transition region in 5 = AB, and 
found AB ~ uj 1 I uz , in agreement with (|4.4j ). Further, the u> dependence of AB saturated for 
uj < ksT/h, in agreement with the ideas we have discussed here |2^J^]. However, Engel et 
al. did not analyze the u and T dependence of a precisely at the critical field 5 = 0. 
It appears to us that it should be relatively straightforward to extend their measurements 
at 5 = to test the validity of the scaling form (|1.2j). Further, it should also be possible to 
theoretically determine the expected form of the scaling function S in the models for the 
quantum Hall transitions considered in Refs 
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A limited test of (|1.2j ) should also be possible in measurements of the a.c. conductivity 
at relatively low frequencies at which Tiuo <C k^T . At such frequencies, we can expand the 
scaling function S(cJ) about To = 0, and analyticity of T > properties then implies the 
form 




(d-2)/z 



S(0) 



V)(0) 



+ ... 



(4.5) 



where S(0) and E^ 2 -*(0) are expected to be positive universal numbers of order unity. Exper- 
imentally, one can test if the frequency- dependent correction to a has the 1/T 2 dependence 
predicted by (|4.5|) . 

There do not appear to be any existing measurements of the a.c. conductivity near 
quantum critical points in other systems. The results ([i~2] ) and ( |4.5|) should apply also to 
the superfluid-insulator transition in thin films |l],|2| , to the quantum transition in the doped 
cuprates || and two-dimensional MOSFETs ||10|| , and to the metal-insulator transition in 



three-dimensional Si:P [[Tl],|T2|]; we hope that experiments will be undertaken, as the results 
will be central to our theoretical understanding of these systems. In particular, if scaling as 
a function of hLo/ksT is observed, it would establish quite conclusively that interactions are 
an essential ingredient in the critical theory. 

Also, we note that our picture suggests a rather interesting non-monotonic u> dependence 
of a in at the metal-insulator transition in d — 3: for small ^cu/fc^T a should decrease with 
increasing to as predicted by Q4.5|) , but for larger hio/hsT it should increase as ~ to 1 ^. It 
would be worthwhile to undertake analytic calculations for disordered, interacting, electronic 
systems to search for this non-monotonic to dependence: there are cases where the critical 
theory is accessible at low orders in the 2 + e expansion (Tl| . There are existing calculations 
for the to and T dependence of the conductivity in the weakly disordered metal |30 



in 



which l/Tip -C fceT/fr; these need to be extended to the critical point where we expect that 
1/% ~ k B T/h. 



E. Self duality 

Of great current interest is the issue of "boson-vortex duality" at two-dimensional quan- 
tum critical points. At the superfluid-insulator transition, many of the experimentally mea- 
sured values of E(0) appear to be tantalizingly close to a value predicted by self-duality 
arguments fl6yi~7fl . For the quantum Hall transitions, experimental evidence for self-duality 
has been presented recently p7J. This self-duality appears rather surprising as there is no 
fundamental reason for an equivalence between the underlying boson and vortex Hamilto- 
nians We suggest here that these inequivalent Hamiltonians will be apparent in the 

value of the T = conductivity, given by £(oo), which, incidentally, is the quantity that 
has been explicitly or implicitly studied in earlier work |T4|--p5|,^2]--|56| . In contrast, the d.c. 
conductivity, given by S(0), is controlled by the crossover to the hydrodynamic, collision 
dominated regime, and we propose that this could be less sensitive to the details of the 
boson-vortex Hamiltonian. This proposal is easily subject to experimental tests: measure 
S(oo) by determining the conductivity in the regime fuo ^> ksT, and see if the self-duality 
predictions continue to hold — we predict they will not. 
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Our suggestion that E(0) could be insensitive to differences in the interaction Hamil- 
tonian of bosons and vortices is motivated by the expectation that an important role of 
interactions is in the collisions that establish local thermodynamic equilibrium p3|j68|| . As 
a result, the equations governing the net hydrodynamic flow of bosons and vortices in the 
collision dominated regime could be more symmetrical than the underlying Hamiltonian. 
In other words, we are proposing here that a true understanding of the experimentally ob- 
served duality near the quantum critical point will emerge from a study of the crossover from 
the microscopic quantum-critical physics of the elementary excitations to the low frequency 
collision dominated regime best described by a quantum Boltzmann equation. 

We suggest that even the simple model S is self-dual for d = 2, n = 2 in its d.c. transport 
i.e. 27rE(0) = 1 exactly for the two dimensional quantum XY rotor model. This model has 
boson particle/hole excitations with short-range interactions, while the vortices in the dual 
representation have long-range logarithmic (Coulomb in d = 2) interactions. However, at 
non-zero T, this logarithmic interaction will be screened in manner analogous to the classical 
Debye screening above the Kosterlitz-Thouless temperature, T KT ; indeed the region g = g c , 
T > of interest here is continuously connected to the g < g c , T > Tkt region. So the 
effective interactions between the bosons and vortices are both short-range for T > 0, leaving 
open the possibility of self-dual behavior at low frequencies. There is now general agreement 
that 27r£(oo) « 0.3 for this model ||19|| . This paper contains the first computation of S(0), 
and as we noted earlier near Eqns ( |1.4| , |1.5|) , it is quite remarkable, though possibly fortuitous, 
that our leading order result for E(0) in the e expansion differs from the self-dual value by 
less than 4%. Definitively establishing this self-duality would however require techniques 
other than expansion in e = 3 — d, or 1/n, as it is only possible precisely at d = 2, n = 2. 
It would be of great interest to undertake higher precision quantum Monte Carlo, exact 
diagonalization, or high temperature series studies while carefully examining the reliability 
of the analytic continuation to real frequencies. 
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Note added: The reader may be interested in a recent article by one of us discussing 
this work in the context of recent results on non-zero temperature dynamical properties of 
S for other values of d and n, including exact solutions in d — 1. 



23 



APPENDIX A: TRANSPORT IN DISSIPATIVE QUANTUM MECHANICS 



We will consider the transport properties of a well-studied and representative model from 



the subject of dissipative quantum mechanics and related quantum-impurity problems |p6 
Our purpose here is to contrast the u and T dependence of the transport coefficient in such 
a situation, with that of bulk critical points studied in the body of the paper. We believe 
that such an exercise will help clarify the significance of our results for the reader. We thank 
E. Fradkin for posing the questions which led to the analysis below. We will use units with 
h = k B = 1 in this appendix. 

Consider the motion of a quantum 'particle' in a periodic potential in the presence of a 
linear coupling to a Ohmic heat bath |36| , |70| -fr^,^8|-^0| . We represent the time-dependent 
co-ordinate of the particle by X(t), and its analytic continuation to imaginary time by X(t), 
and the Fourier transform to imaginary Matsubara frequencies by X(uj n ). The imaginary 
time effective action obtained after integrating out the degrees of freedom of the heat bath 
is 

S i = Kl \ X Mf -yj drcos(X(r)). (Al) 

Here a is a dimensionless coupling constant characterizing the strength of the Ohmic dis- 
sipation, and y measures the strength of the periodic potential. Models like ( |A1| ) describe 
tunneling in a SQUID, where X is interpreted as the flux in the SQUID, or tunneling be- 
tween Luttinger liquids and quantum Hall edge states [^.^1].^,^], where X now becomes 
a bosonic phase field. 

The action ( |A1[ ) can be written in a form local in time if we extend the quantum degree 
of freedom X(t) to an infinite number of degrees of freedom X(x,t) lying along the line 
-oo < x < oo, with X{t) = X(x = 0,t); then Si is equivalent to |7l],|f|,|7|,|76|,|38|-|g] : 



ac 
8^ 



J dxdr ((d x X(x,r)) 2 + ^(d T X(x,r)) 2 ^ - y J rfrcos(X(r)), (A2) 



where c is an arbitrary velocity. Notice that the cosine interaction acts only along the single 
line x = 0, identifying this as a boundary critical phenomena problem. 

We shall be interested here in the linear response of the system to a time- dependent 
force, F(t), acting on the particle. In imaginary time, in the presence of such a force 

Si -> Si - J cItF{t)X{t). (A3) 

In the presence of a time-independent force, we expect the particle to acquire a finite velocity, 
V = dX/dt, in steady state: this allows us to define a mobility, G, by V = GF. More 
generally, we expect a frequency dependent response G(co), defined by 

V(u) = G{uj)F{u). (A4) 

It is our purpose here to describe the behavior of the dynamic mobility, G(u), at low u and 
T, and to compare it with the results for <j(lv) obtained in the main part of the paper. We 
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also note that G(u) is the conductance associated with tunneling between quantum Hall 
edge states ||(| at a value g = 1/a in the notation of Moon et al. [fnj. 

Consider first the properties of S\ under a T = renormalization group transformation 
under which r — > re~^. It is found |fF0H72]] that a remains unrenormalized, while the 



potential strength y obeys the simple flow equation 



dy 
di 



a 



-is 



(A5) 



For a ^ 1, this flow has fixed points only at y = and y — oo (we will not consider the 
case a = 1 here). For a < 1, the y = fixed point is stable and the y = oo fixed point is 
unstable; for a > 1 the stability of the fixed points is interchanged. Notice that both fixed 
points are free field theories, supplemented by free or fixed boundary conditions at x = 0. 
This is a crucial difference from the bulk theory S in which the critical theory is interacting 
for d < 3, and is primarily responsible for the differences in the structure of G(u) and o~{u) 
we shall find below. 

Let us now write down the general scaling predictions for G(u) which follow from the 
renormalization group arguments. The system will be completely characterized by a single 
nonuniversal energy scale, Tk, which measures its deviation from the unstable fixed point. 
From the flow equations (|A5|) we can deduce Tk ~ for a > 1; the perturbation 

theory to be discussed below shows Tk ~ also for a < 1. The mobility G has a 

zero scaling dimension, and therefore obeys the scaling form 



G(oj) 



uj T 



(A6) 



where Q a is a universal function. We now describe the low T, uj form of Q a for the cases 
a < 1 and a > 1 separately. 

Consider first a < 1. In this case y flows to zero, and therefore the periodic potential has 
vanishing strength at long times and the particle is delocalized. Naive perturbation theory 
in y is expected to be reliable. The mobility can computed in this perturbation theory using 
a Kubo formula: such a computation was carried out in Appendix A of Ref J75|, and we 
find from their results 



r>( \ 271 ( T 
G(u) = — 

a \1k 



2(l/«-l) 



where the ellipses represent terms which are higher order in T /Tk — this result 
uj <^ T K , T <^ T K , but u /T can be arbitrary. Notice that ( |A^ ) is consistent with ( 
explicit form of V a can be deduced after some analysis of Eqn (A7) of Ref ]7fJ : 



(A7) 

holds for 
. The 



VJuJ) 



, 1 iuj 

a Z7T 



sinh(a;/2) 

UJ 

tan(7r/a) 
ZJ 



r 



1 iuj \ 
a + 27r/ 



cosh(u;/2) 
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/n s 2(l/a-l) 

V a (uJ ^ oo) = (-) 2 1 - 2 / a [l-ztan(7r/a)] 



7T, 
1 cJ 2 

*Pi/2(p) = g + ^""2 ' a s i m pl e special case |75| , |76| . (A8) 

The overall normalization of V a is arbitrary, as it can be absorbed into a redefinition of Tk- 
The most important property of (|A7|JA8|) is that the leading term in the mobility, 2n/a, is 
independent of uj/T. Alternatively stated, we have 

2-7T 

G(u -> 0, T = 0) = G(uj = 0, T -f 0) = — . (A9) 

However, it is clear that this leading term is a property of the y = fixed point, and the 
independence on uj/T is a consequence of its free field nature. The next, subdominant, 
term arises from the leading irrelevant operator in the theory, and does indeed lead to a 
non-trivial dependence on uj/T in the universal function V a . In contrast, the bulk model S, 
not being a free-field at the critical point, had such a uj/T dependence already in the leading 
term, without the inclusion of any irrelevant operators. This distinction is the central point 
of this appendix. 

Now turn to a > 1. In this case y flows to oo, implying that the periodic potential 
localizes the particle at T = 0. The low T mobility clearly cannot be computed by a naive 
perturbation theory in y. Instead one can use a self-duality property of Si |70| , [72| , |76]| under 
which a «-> 1/az, and then use perturbation theory. The y = oo fixed point now implies that 
the leading scaling result for the mobility is simply G = 0, and so the uj — » and T — > 
limits commute again for a trivial reason. The leading uj and T dependence is given by an 
irrelevant operator, which now yields 

Gip) = (jr) Pi/a (£)+■■■> ( A1 °) 

where the scaling function V\/ a was defined in (|A8|), the corrections are higher order in 
T/Tk, and the result holds for uj <^ Tk, T Tk, but uj/T arbitrary. Notice also the result 

G(oo -^0,T = 0) = G(lu = 0,T ^0) = 0, (All) 

which arises from the y = oo free field fixed point. 

We also note for completeness the high T behavior, with T ^> Tk, to allow contact 
with the quantum Hall edge state tunneling results of Ref J77J. The key is to note that 
the a <-> 1/a duality interchanges small and large y. So the result ( |A10|) also describes the 
T ^> Tk limit of the a < 1 case, while (|A7|) describes T 3> Tk for a > 1. The main change is 
that the values of Tk (denoted Tk) inserted in these expressions will not be the same as those 
in the low T limit. Perturbation theory cannot determine the universal relationship between 
Tk and Tk- this requires use of the exact integr ability of Si, as discussed in Refs [j^,|73|-|75|l . 

Finally, we suggest that the above commutativity of the limits uo — >• and T — > (and 
also uj — > oo, T — > oo, where oo refers to a scale much larger than T K ), and the free field 
nature of the fixed points may be the reason for the success of the imaginary time Monte 



Carlo simulation of Moon et al. |77]] 
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APPENDIX B: GENERAL TRANSPORT EQUATION 



Here, we will generalize the transport equation (|3.25|) of the model S to arbitrary n > 2 



First consider 0(n = 3). This case applies to quantum antiferromagnets |H| and the 
external 'potentials' U a (a = 1, 2, 3) correspond to the three components of a magnetic field. 
Let us take the field pointing in the 3 direction: U a = (0,0, H). For the antisymmetric 
matrices we choose L^g = € aa/ 3, the third-rank antisymmetric tensor. Note that such a field 
is the same as that in ( |2.1| ). The distribution functions become diagonal by changing basis 
from £11,2,3 to 

MM S ai(M) ^ aa(M) .MM) (Bi) 

The particle distribution function now has the three diagonal components f±, / 3 . The 
'magnetization' current is non-zero only in the a = 3 direction and equals 

J3 = Q I 7^r[/+(M-/-(M)] 

J (2-n) a e k 

A key simplifying feature is that the distribution function f 3 must be even in the external 
field, and will therefore only be modified at quadratic order in H. We will be satisfied by 
working in linear response, in which case / 3 (k, t) = n(e k ). 

The n > 3 case is very similar: the only difference is that there are now n — 2 values of 
a for which / a (k, t) = n{ek) in linear response. 

The generalization of the transport equation ( |3.25| ) to 0(n) is now easily obtained by an 
application of Fermi's golden rule: 

(9_ xn ^ R d\ f(v ,,_uh[ d% d% d% 1 

\dt + ^ dk J JX[ ' ' ~ 9 J (2ny (2nY (2ttY 16e k e kl e k . 2 e k3 

x (27r) d 5(k + ki - k 2 - k 3 )2n5(e k + e kl - e k , 2 - e k3 ){ 

4/ A (k, t)/_ A (k!, t) [1 + / A (k 2 , t)] [1 + /_ A (k 3 , t)} 

+ 2/ A (k, t)/ A (k ls t) [1 + / A (k 2 , t)\ [1 + / A (k 3 , *)] 

+ {n- 2)/ A (k, t)n{e kl )[l + / A (k 2 , t)][l + n(e fc3 )] 

+ ^ y ^/A(k,t)/_ A (k 1 ,t)[l + n( £fc2 )][l + n( £fc3 )] 

- 4[1 + / A (k, t)] [1 + /_ A (kx, t)]/ A (k 2 , t)/_ A (k 3 , t) 

- 2[1 + A(k, t)] [1 + / A (k l5 t)]/ A (k 2 , t)/ A (ka, t) 

- (n - 2)[1 + / A (k,t)][l + n( £fel )]/ A (k 2 ,t)n(£ fe3 ) 
fn-21 



■[1 + / A (k,t)][l + /_ A (k 1 ,t)]n(£ fc2 )n(£ fe3 )} (B3) 



with A = ±1. The analysis of the linearized form of this equations is similar to that for 
( |3.25| ), but will not be presented here. Such an analysis will lead to a determination of 
the spin conductivity, a, which is related to the spin diffusion constant D s by the Einstein 
relation 
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o = D sX, 



(B4) 



where x is the uniform spin susceptibility. Results for the x a t non-zero T above the critical 
point have been given earlier in the 1/n pT| , |i"3"f and e expansions. 



APPENDIX C: COMPUTATIONS WITH THE COLLISION TERM 

We describe here some of the steps in between the original transport equation Q3.25| ) and 
the linearized form (|3.28|) 



First we insert the parametrization ( 3.22 ) into ( p.25|) , and linearize in the the external 
electric field. Then, notice that in the collision term the unknown function ip appears only 
in the integrals over the radial components of the momenta. The angular integrals involve 
only known functions and can be performed analytically. As already noted in Section |111B 



the integrals in the collision term can be done directly in d = 3 to obtain the leading result, 
and so all the computations here are for d = 3. One needed angular integral is 



J dk 3 g(k 3 ) J dtt kl dtt k2 dtt k3 5 3 (k + k 3 -k 1 - k 2 )8(k + k 3 — k 1 - k 2 ) 

-gih + h-^hi^kuh), (CI) 



8tt 2 



kk\k 2 {k\ + k 2 — k)' 

where g(k 3 ) is some function of k 3 . The result involves the function I\ which is given by 



— k k\ + k 2 > k and k\ < k and k 2 < k 

(C2) 



h(k, h,k 2 



k\ + k 2 < k 

k\ + k 2 — k k\ + k 2 > k and k\ < k and k 2 < k 

k\ for k\<k and k 2 > k 

k k\ > k and k 2 > k 

k 2 k\ > k and k 2 < k 



A second angular integral has a single vector momentum in the integrand 
J dk 3 g(k 3 ) J dfi kl dfi k2 dfi k3 k 1( 5 3 (k + k x - k 2 - k 3 )5(k + k x - k 2 - k 3 ) 

87T 2 k 



3fc 3 fcifc 2 (fc + h -k 2 ) 



gik + h-k^hi^kuh), (C3) 



where now 



h(k,ki,k 2 ) = < 



k 2 >k + ki 

k 3 + k\ + 2kl - 3kl(k + h) + 3kk x k 2 k 2 <k + k x and k 2 > k and k 2 > k x 

k 3 for k 2 > k and k 2 < ki 

—2k 3 + 3k 2 (k + ki) — 3kkik 2 k 2 < k and k 2 < k\ 

kf k 2 > k\ and k 2 < k 

(C4) 



Finally, we also need a variation of (|C3|): 
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J dk 3 g(k 3 ) J rffi kl rffi k2 dfi k3 ki5 3 (k + k 2 - k x - k 3 )5(k + k 2 - h - k 

8vr 2 k 



3k 3 kik 2 (k + k 2 — ki 



g(k + k 2 -k 1 )h(k,k 1 ,k 2 ), (C5) 



where 



I 3 (k, k x ,k 2 ) 



Using 
becomes 



k x >k + k 2 

k 3 — kf — 2k 2 — 3k 2 (k — ki) + 3kk x k 2 k\ < k + k 2 and ki>k and k 2 < ki 

k 3 for k\ > k and k 2 > k\ 

-2k\ - ?>k\{k -h) + 3khk 2 k x <k and k 2 < k x 

k\ k 2 > k\ and k\<k 

(C6) 

in the linearized form of fl3.25p , the integral equation satisfied by ip 



. . ldn(k) 

—t<jip(k, on) + — 



k dk 
tie 2 \ 18^j(k,uj)k 2 roc 



75fc 4 



n(k) 



o 



dkidk 2 Ii(k, ki, k 2 )n{k 2 )n{ki) [1 + n{ki + k 2 — k)] 



+ 2[1 + n{k)} I dk x dk 2 ^ k ^ I 2 {k, h, k 2 )n{k 2 )n{k + k x - k 2 ) 

-An{k) I dk x dk j^ Y- h{k, h, k 2 )n{k 2 )[l + n{k + k 2 - fci)]| . 

nyki) J 



(C7) 



This can be turned into the form (|3.28| ) by evaluating the integrals over k 2 , which can 
also be done analytically To do this, it is first necessary to separate the products of two 
Bose functions each with k 2 in their arguments, into terms in which only one Bose function 
involves k 2 . this is done by repeated use of the identity 



n{k 2 + a)n(k 2 + b) = n(b — a)n(k 2 + a) + n(a — b)n(k 2 + b) 



(C8) 



which is valid for any a, b. From the forms of the functions h,2,3 above, it is now clear 
that we only need integrals of k 2 over Bose functions of k 2 , times powers of k 2 . The general 
integral needed is of the form 



dko 



An 



; (fc 2 +c)/T _ l 



dk? 



h 2 



o{k 2 +c)/T 



dk- 



h 2 



e (k 2 +c)/T _ I 



(C9) 



where n is an integer, and c is a real number. Changing variables on the right hand side, 
we can rewrite the first term on the right hand side as 



dh 



7,71 
"2 



e (fc2+c)/T 







00 d u (^ + Q ) ? 

aK2 e (k 2 +c+a)/T 



(CIO) 



and similarly for the second term. Expanding the polynomial in the numerator, we finally 
conclude that all of the integrals over k 2 can be reduced to the following basic integral: 



29 



l dk \( k2+ a)/T_ l = Tn+l nn + l)Li n+1 (e-^). (Cll) 
Here Li p (z) is the polylogarithm function of order p, defined by the series 

Li P W = E^- ( C12 ) 

n=l '* 

Notice Lip(l) = ((p). 

After performing the integrals over hi as described above, (|C7|) takes the form ( |3.28| ). 



We will not display explicit expressions for Fi :2 (k, k\) here as they are quite lengthy and not 
particularly informative. We used a CERNLIB routine for numerical evaluation of IA p (z) 
for < z < 1: this allowed very rapid determination of the kernel of the integral equation 

(ED- 



APPENDIX D: ORDER PARAMETER SUSCEPTIBILITY 



This appendix will review the results of Ref. |29| on the finite temperature crossovers in 
the order parameter susceptibility \ — ((fia'fia) (no summation over a) 



Ref. 



was con- 
cerned with making statements correct to all orders in e and to understanding the crossover 
to classical critical fluctuations in the vicinity of the finite temperature phase transition 
line, and this required a rather involved analysis. Here, we are mainly interested in low 
order results at finite temperatures above the quantum-critical point and the quantum dis- 
ordered phase: in this case the necessary results can be obtained directly in a self-consistent 
Hartree-Fock-like analysis, as we now show. 

The self-consistent Hartree-Fock susceptibility of S is obtained by summing all the one- 
loop tadpole diagrams; this leads to the expression 



u 2 n + k 2 + m 2 (T) 
where m 2 (T) is given by the solution of the equation 



PI) 



m 2 (T) = t + < + Mo (!^)rE/ 



d d q 



6 



(2vr )^2 +g 2 + m 2 (T) 

d d q 1 r d d+] 



(27r) d e 2 n + q 2 + m 2 {T) 



p 1 



(2 



7T 



(D2) 



In the second equation we have inserted the leading result for the value of mg c . We are 
following the convention here of denoting d + 1 dimensional spacetime momenta by p, and 
d dimensional spatial momenta by k, q. In the critical region, the coupling uq is of order e, 
so one might think that it is permissible to set m 2 (T) = t on the right hand side of ( |L)2| ) 
and obtain a result correct to order e 2 . However this is not adequate for our purposes for 
two reasons [29]: 

(i) The resulting expression for m 2 (T) is not analytic as a function of to at to — for T > 0. 
This analyticity is required on the physical ground that there can be no thermodynamic 
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singularity at T > at the quantum-critical coupling t = 0. 

(ii) At T > above the quantum critical coupling t there is a contribution to m 2 (T) of 
order e 3//2 which is missed. 



To proceed, we can either use the analysis of Ref , or directly analyze the singularity 
structure of ( [D2j ) as u — ► 0. By either method, it can be shown that it is permissible to set 
m 2 (T) = to in only the e n 7^ terms on the right hand side of (|D2| ). To describe this, we 
write 



m 2 (T) = R{T)+5m 2 {T) 



where 



R(T) =t + u 



n + 2 



d d q 
(2tt) c 



+ 



T 



d d+1 p 1 
(2vr) a! +V 



(D3) 



(D4) 



was a quantity introduced in Ref p9|, and 



5m 2 (T) = u ( 



n + 2 
6 



d d q ( 1 1 
(2vr) d l v g 2 + m 2 (T) ~ ^ 



(D5) 



Finally, on the right hand side of (|05|) we replace m 2 (T) by R(T). The resulting expression 



for Sm (T) when inserted along with (|D4j ) into ( P3|) gives us our final result for m (T). 

We will now manipulate the above result for m 2 (T) into a form which explicitly displays 
its universal scaling nature. First replace the bare couplings by renormalized couplings; to 
the order we are computing things, this is equivalent to performing the substitutions |79|j80| 
Mo = (J, e g I 'Sd+i and to = t(l + (n + 2)g/(6e)) where /x is a renormalization momentum 
scale, g is a dimensionless coupling constant, and t is a renormalized coupling measuring the 
deviation from the critical point. Expanding the result to order g, we find that the poles in 
e cancel. Finally, we set g at its fixed point value [ff9|JS0| g = g* = 6e/(n + 8). Our final 
result for m 2 (T) then becomes 



m 2 (T) = R(T) 



n + 2 
n + 8 



2ttTJR(T) 



(D6) 



where the second term on the right hand side is the contribution of Sm 2 (T). A somewhat 



subtle analysis, discussed at some length in Ref |[29|| , is required to obtain the following final 
form for R(T): 



i?(T)=t(l + e— ln- 
with the crossover function G(y) given by 



+ 2 

+ eT 2 



n + 8° \T- 



t 



(D7) 



G{y) = -2 dq 



In 2q 



2V cosh(vV + y) 



q 2 + y 



q 



2Jq 2 + 1/e 



;ds) 



This form of G(y) is valid for both negative and positive y (when the argument of the square 
root is negative we use the identity coshix = cosx) and is easily shown to be analytic at 
y = where 
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2tt 2 

G(0) = , dG/dy(0) = 2.453808582 .... (D9) 

3 

The result ( p7[ ) can therefore be used on both sides of the critical coupling t — 0, and the 
required analyticity of T > properties at the critical coupling has been achieved. Also 
notice that at t = 0, R(T) is of order e, so the result ( p6| ) contains a term of order e 3 / 2 . 

Finally, let us express m 2 (T) in terms of experimentally measurable energy scales. The 
energy scales have to be defined differently for t > and t < 0. 

For t > we choose the T = energy gap, A, to measure the deviation from the critical 
point. Using the relation |29f A 2 = p 2 {t/ p 2 ) 2u in flD7|) (i/ = 1/2 + e(n + 2)/(4(n + 8)) is the 



correlation length exponent) we get our final universal expression for R(T), valid everywhere 
above the quantum-disordered (insulating) phase (t > 0) 



-R(T) = A 1 + e — ' — In — + eT — ' — G (D10) 
1 n + 8 A/ Ti + 8 1 ' 




Notice that the arbitrary momentum scale /x has disappeared, and combined with ( |D6| ) we 
now have a result for m 2 (T) expressed solely in terms of the measurable energy scales A 
and T. 

For t < 0, our results are confined to the normal phase T > T c (t). We measure the 
deviation from the critical point by the value of the superfluid stiffness, p s at T = (for 
n = 3 p s is the spin stiffness of the ordered antiferromagnetic phase). To obtain a quantity 
with the dimensions of energy we define 

p s = ((3 - d) Ps ) l/{d - X) . (Dll) 

The numerical factors are for future convenience; also notice that in d = 2, p~ s = p s . Then 



we use the expression for p s in terms of t and p in Ref |29] to express t in terms of p s and 



p in (p37|) . Then all the p dependencies cancel as before, and we get our final result 

~2 

R(T) = P»- I 1 



n + 8\ 2(n + 




i2 n + 2 



+ eT — ' — G —-. — ^—^\ (D12) 

Again combined with (|D^) we now have an expression for m 2 (T) in terms of the superfluid 
stiffness of the ordered phase at T = for t < and T > T c (t). 
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1 — ► 

CO 

FIG. 1. A sketch of the expected form of the real part, of the universal scaling function £ 
appearing in the scaling form (|1.2| ) for the conductivity, as a function of To = huj/ksT. There is 
a Drude-like peak from the inelastic scatterings between thermally excited carriers at u; of order 
unity. At larger ZJ, there is a crossover to the collisionless regime where £' ~ zj( d ~ 2 )/ z as uJ — > oo. 
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1 ► 

CO 

FIG. 2. Universal form of the conductivity <r(u, T — > 0) in d = 2; the vertical scale is measured 
in units of h/Q 2 . Only the uj = value is given by the universal number S(0). For all cj > 0, 
(h/Q 2 )a = £(oo). Q is the 'charge' of the order parameter: for the superfluid-insulator transition 
Q = 2e, while for quantum antiferromagnets Q = g\xn- 
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FIG. 3. Structure of the real part, of the universal scaling function £ in ( |1.2| ) for the conduc- 
tivity at the quantum critical coupling of the model S defined in ( |1,3[ ) . The spatial dimensionality 
d = 3 — e, and e is assumed to be small. As before ZJ = fiw/kBT. The Drude peak at small ZJ has 
a width of order e 2 and a height of order 1 /e 2 : this feature of the conductivity is denoted later in 
the paper by 07. The collisionless contribution (denoted an later) begins at ZJ of order e 1 / 2 ; as 
ZJ — > 00, this contribution is a number of order unity times ZJ 1 ^ 6 
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FIG. 4. The real part, S', of the universal scaling function E (see (|1.2|) ) for the conductivity 
at the quantum critical coupling of the model S, correct to first order in e = 3 — d. The numerical 
values are obtained from and ( |2~To| ) with d = 2 (e = 1). There is a delta function precisely 
at uj/T = represented by the heavy arrow: the weight of this delta function is given in ( |2.7D and 



(2.11). The delta function contributes to 07, and the higher frequency continuum to an 
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co=o 




FIG. 5. 



Real part of the universal function k^^{k,u>) as a function of k for a few values of 



. The function is defined in ( |3.30| ) and ( 3.22| ), and was obtained by numerical solution of the 
linearized quantum Boltzmann equation fl3.3ip . At u> = 0, ^ is real, but is complex for general ui. 



Here k = k/T, and ui = uj/e 2 = Lo/e 2 T (in physical units k = hck/kgT, lo = hu>/e 2 kBT). 
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FIG. 6. The real part of the universal function S/(S), which is related to the low frequency 
part of the conductivity (aj(oo)) by Q3.32 ). The results are obtained by the numerical solution of 
( 3.31| ), followed by the integration in fl3.33| ). This function describes the inelastic collision- induced 
broadening of the uj = delta function in Fig ^| at a frequency scale of order e 2 T. The conductivity 
has an additional continuum contribution (au(uj)) at frequencies larger than uj ~ e x / 2 T which is 
not shown above (see Fig H). 
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